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Abstract 



We study the influence of disorder on propagation of waves in one-dimensional struc- 
tures. Transmission properties of the Schrodinger equation with the white noise potential 
can be expressed through the Lyapunov exponent 7 which we determine explicitly as a 
function of the noise intensity a and the frequency to. We find uniform two-parameter 
asymptotic expressions for 7 which allow us to evaluate 7 for different relations between 
a and to. 

> 

in 1 Introduction 

m 



Analysis of the Lyapunov exponent plays an important role in the problems of propagation 
of waves. Much attention in recent years has been given to periodic and disordered one- 
dimensional structures (see, e.g., p] and references therein). Application of such structures 
include photonic crystals, coupled-resonator optical waveguides, slow-light devices, metamate- 
rials, semiconductor superlattices with given transmission properties and more. Characteristics 
^ \ of real materials, however, may be considerably different from predicted ones due to inevitable 
disorder of geometric and material parameters. Structural perturbation and manufacturing 
inaccuracies may lead to a substantial reduction of performance [2] , [3] and even to complete 
suppression of wave propagation [I], [5]. In the strongest form the disorder may manifest itself 
in the Anderson type localization, which is a general wave phenomenon for classical and quan- 
tum waves in disordered systems. Due to multiple scattering, destructive interference occurs, 
so that the wave amplitude decreases exponentially in the medium. The rate of decay of the 
solution is defined by the Lyapunov exponent, 7 (the inverse of the localization length). 

The Lyapunov exponent is of major importance in the problems of slowing down light by 
orders of magnitude, a phenomenon which has been extensively discussed in the literature [6]. 
Such a possibility can be useful in a variety of optical and microwave applications. Choosing 
wave frequency near stationary points of the electromagnetic dispersion relation allows one to 
reduce the group velocity by two orders of magnitude [6]- [7]. Another model for slowing down 
light proposed in [8] is based on the periodic necklace-type waveguide. An essential common 
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feature of such devices is that the frequency of propagating wave must be located near the 
spectral band edge or degenerate band edge of the corresponding medium. 

One of the major factor, however, limiting performance of optical waveguides is their struc- 
tural imperfection, and investigation of the influence of different kinds of disorder on the Lya- 
punov exponent becomes of primary importance in manufacturing of optical devices. Papers 
[H]-[in] performed sensitivity analysis of photonic crystals near degenerate and regular spectral 
band edge under perturbations of geometrical and material parameters, and the impact of these 
perturbations was compared. Numerical simulations were used to investigate the sensitivity of 
field enhancement effects at Fabry-Perot resonances against the bandwidth of the excitation. 
The gigantic amplitude increase produced by such resonances was observed, which dramati- 
cally deteriorates under small perturbations of the thickness of layers. In particular, the peak 
field intensity of a 16 layer stack regular band edge photonic crystal decreases by an order of 
magnitude under 1% perturbation of the layer thicknesses. 

The central idea of chaotic behaviour of dynamical systems is widely seen as strong de- 
pendence on initial conditions [UJ. In a chaotic system, a very small change in the initial 
conditions can cause a big change in the trajectory. Quantitative criterion for sensitivity is 
characterized by the Lyapunov exponent (also called the maximal Lyapunov exponent) which 
measures exponential growth rate of the solution. 

In all the applications, the influence of the strength of perturbation o is essentially different 
depending on how close the frequency uj is to the band edge. While 7 ~ o 2 in the bulk of 
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the band, it has order 7 ~ as near the band edges. Therefore, it is important to construct 
uniform two-parameter asymptotic expansion of the Lyapunov exponent for all frequencies of 
wave propagation. 

Let 7j = 7x(wi, u) be the Lyapunov exponent of the optical equation 

— ip"(x) = uj\n 2 (x, a)xp(x), (1.1) 

where < a ^ o~o denotes the strength of perturbation of either the value of the periodic 
refraction index n(x, 0) or the length of its period, and uj± is located in an ^-neighbourhood of 
a non-degenerate band edge. Let also 7 = 7(0;, a) be the Lyapunov exponent of the Schrodinger 
equation 

— if)"(x) + aw(x)%jj(x) = uj 2 ip(x), (1.2) 

with the white noise potential w(x) and the frequency u located near the spectral edge u = 0. 
The frequency uj is defined by u\ using the relation 

cos uj = ^ tr M (coi), Ui G Q. (1.3) 

Here M is the transfer matrix of ( 11. ip . Then in [12] it was shown that r y 1 and 7 estimate each 
other uniformly. That means 

C17 < 7i < c 2 7, (1.4) 

where positive constants c\ and C2 depend only on cr an d ^ but not on the values of a and 
uo\. Thus, a uniform two-parameter asymptotic expansion of the Lyapunov exponent of the 
Schrodinger equation (II. 2p provides a uniform estimate of the Lyapunov exponent of the wave 
equation (11.11) . Because of this, we consider in details the Schrodinger equation with different 
relations between the intensity of noise a and the frequency uj. We will obtain a uniform 
asymptotic expansion of 7 followed by numerical evaluation of 7 to determine the range of 
parameters for which the asymptotic formulas can be used. 
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2 Integral representation of the Lyapunov exponent 



In this paper we study the asymptotic behaviour of the Lyapunov exponent 7(A, a) of the 
random one-dimensional Schrodinger operators 

Hip = —if>"(x) + aw(x)ip = \ip, x > 0, 
ip(0) = cos 6*o, = sin 6*o, 

CO 

as a function of the energy A and the coupling constant (intensity) a of the white noise potential. 
Here uj = y\X\ (equations with any real A will be considered). The white noise w(x) can be 
understood as the limit of a stationary Gaussian process with independent increments on 
intervals of length A, A < 1 (cf. (12]): 

w(x) = hm ~~^j=^i x e [0) A)> and is Af(0, 1) random variable. (2.2) 

An equivalent alternative approach to define the operator with a white noise potential (which 
is used in the present paper) is based on the reduction of (12. ip to a system of Ito stochastic 
differential equations of the first order for the phase 9 and amplitude r of if). In what follows 
we will denote z = — cot 9. Then the Ito and Stratonovich representations of (12. ip coincide in 
coordinates (z,r) and both give the same result (see [IB]). 
After the standard substitution (in the Priiffer form) 

u = ip, v = —, cu = V A , 

equation (12. ip takes the form 

du = ljv dx, 

, audw (2.3) 
dv = —uju dx H . 

UJ 

Usually the reduction of the Schrodinger equation to a system of two first order equations is 
accompanied by the phase-amplitude formalism, i.e., by using the polar coordinates 

u — if) — r(x) sin 9(x), 

iP' (2-4) 
v — — = r(x) cos9(x), 
cj 

where 9 G [0, 7r). 

The Lyapunov exponent 7(A, a) is defined by 

. . , lnr(x) 

7 A, a = lim — ■ 2.5 

t-yoc x 

This limit exists P-a.s. on the probability space (ft, J 7 , V) defined by the white noise w(-) 
in (12. ip (see [13], [E]). We will discuss the limit behaviour of j(\,a) in the different regimes 
when a — > or a — > oo (small and large disorder) and = u — > or yfAj = u — > oo (long 
and short waves). We also introduce some combinations of A and a to describe the intermediate 
asymptotics. 
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Asymptotic behaviour of 7 for a — > and fixed A > was established in papers [15J and 
[16] where it was shown that 

7~|^, ^0. (2.6) 

Paper [T7] contains the formula 7(0, a) ~ as for fixed A = and a — > 00. We show that in fact 
7(0;, cr) = usf(v), where v = and / is a function of v. Thus, (I2.6P is also valid for arbitrary 
u and u such that v — > 00. 

In this paper we suggest a simple approach to obtaining the asymptotic expansion for 7 as 
v —7- 00 as well as v — > 0. Therefore, different combinations of large and small perturbations 

2 

and frequencies are covered. We show that 7 = CV3, v — > 0. In particular, it provides low 
frequency asymptotics of 7 when a is fixed and A — > 0. The case of A < is also considered. 
The next section contains an exact formula for 7 followed by a section with derivation of its 
asymptotics expansion. We will use equations (12. 3p to obtain the integral representations for 

7- 

Theorem 1. The following formulas for the Lyapunov exponent hold P-a.s. 
1. IfX>0 then 



where 



1 (u 1 a)=uu- 1 J p(x) ^ 1 + x2 y dx, (2.7) 
p( x ) = Ce*& J e"* w dt, $(x) = v I x + yj , v = -j^, (2.8) 



with a normalizing constant C given by 



fKZ roc -2*(a;) 



2. If A < then 



where 



and 



f°° 1 — x 2 f°° x 

7(w, v) = ujv' 1 I p(x)— — —dx-2u p(x)— — - dx. (2.10) 
J-oc {l + x 2 y 1 + x 2 

p(x) = Ce^ f\-*®dt, = v - x^j , v=^, (2.H) 

[WZ poo -29 (x) 

c-'H-L V fc (2l2) 



Proof. Let A > 0. We rewrite (I2.3P in polar coordinates (12. 4p but use z = — cot 8 instead of 9. 
The resulting equations will be much simpler in the coordinates (z,r) than in (r,8). One only 
needs to keep in mind that the process z(x) = — cot 8(x) has boundary conditions at infinity: 
if z(x) — > 00 when x approaches from the left to some finite value Xq, then z(xq + 0) = —00. 
This jump corresponds to the transition of the angle variable 9 = 6(x) from tc to zero. We will 
use the Ito formula in order to write equations for (z,r). The Ito formula says that for any 
smooth function F = F(u,v), 



dF OF ld 2 F 
dF(u(x), v(x)) = —du + — dv + -—rdv 2 , (2.13) 
ou ov 2 ov 
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where du, dv are given by (12.31) and dw 2 = dx. There is only one quadratic in du, dv term on 
the right, since equation ( 12. 3 p for du does not contain dw. If we put du, dv from (12. 3p into the 
Ito formula we arrive at 

, dF , dF faudw \ 13 2 F(xV , 

dF(u, v) = -r—covdx + coudx ] + --^-^ Trdx 

du dv \ to J 2 dv 2 co 2 

OF OF ld 2 Fa 2 u 2 \ a dF , . 

— ujv - — UU+-— — — —\dxA —udw. (2.14) 

du dv I dv z co z J co dv 

V V 

Since z = — cot 9 = — (see (I2.4p ). we apply (I2.14p to function F(u, v) = z = — . For this 

u u 

function F we have 

dF _ v dF _ 1 d 2 F _ Q 

du u 2 ' dv u dv 2 

and therefore, 

dz = w(l + z 2 ) dx - - dw. (2.16) 

CO 

We will need function In r = | ln(w 2 + v 2 ), not r. Thus from (I2.14p with -F(u, v) = | ln(w 2 + 
f 2 ) we obtain 

O" 2 1 — Z 2 (7 Z 

dlnr = — — — dx -dw. (2.17) 

2co 2 (1 + z 2 ) 2 col + z 2 K J 

The right-hand side of the latter equation does not depend on r. Thus the solution of the 
original Schrodinger equation (12. ip is reduced to a solution of the stochastic equation ( I2.16P 
followed by a direct integration of ( I2.17p . 

Consider the generator of the diffusion process z(x): 

The probability density p(x,Zi,Z2) (with respect to Z2 with initial point z-i) of the diffusion 
process z(x) satisfies the equation 

dp 

— = L zl p = L* Z2 p with p\ x=0 = 8{z x - z 2 ), (2.18) 

and the limiting probability density p is given by 

L*p = Q (2.19) 

or 

The solution of the last equation has the form (12.81) where the normalizing constant is defined 
by the condition 

/oo 
p(z)dz = 1, 
-oo 

that is, 

C- 1 = / e^dz / e-^dr. 



We make the substitution r = s + z here and evaluate the integral in z. This leads to ( 12. 9p . 
After the limiting probability density p(z) is found, formula ( 12. 7p is an immediate consequence 
of (12. 5p and (12.171) . One needs only to note that the boundedness of the integrand in the second 
term in (12.171) implies that P-a.s. 



: dw — > as i->oo. 



If A ^ then we denote (J 1 



l + z 2 

A in (I2.ip . and the system (12 .3p becomes 
du = uv dx, 



dv = ujudx + 



o~u dw 



UJ 



The stochastic differential equation for z(x) = — cot 9 analogous to (I2.16P has the form 



a 



dz = uj(z — 1) dx dw 



id 



The stochastic equation for lnr similar to (I2.17P is 



dlnr 



1-z 2 



2loz 



2lo 2 (1 + z 2 ) 2 1 + z 2 



dx 



a z 



ujI + z 2 



dw. 



(2.21) 



(2.22) 



(2.23) 



Determining the limiting probability density p as in (12.201) we obtain that it satisfies the equation 

(2.24) 



a 2 d 2 p d , 2 v 



2u 3 dz 2 dz 

The solution of this equation is given by (12.111) with a normalizing constant determined from 
(I2.12p . Therefore, the expression of the Lyapunov exponent is given by (12.101) . □ 



3 Asymptotics of the Lyapunov exponent 

The exact expression of the Lyapunov exponent 7 is given by (12. 7p if A ^ and by (I2.10p if 
A ^ 0, and can be calculated numerically. Although equation (12. ip depends on two parameters 
A and a, formula (12.71) implies that 7 depends on the product of u = a/|A[ and a function of 
v = The next theorem provides the asymptotic behaviour of the Lyapunov exponent for 
v — )• and v — > 00. 

Theorem 2. Lyapunov exponent of the Schrodinger operator k2. 1\) has the following asymptotic 
behaviour: 

(a) If X jS; and v — > (long wave asymptotics) then 

7 = c<r§ (l + 0(z/§)) , (3.1) 

where exact value of c is given by Ii3.13\) . 

(b) If X > and v — > 00 (short wave asymptotics) then 
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(c) If A < and v — > oo then 

7 = w(l + 0(^ 1 )). (3.3) 

Proof, (a) Let us find the asymptotic expansion of C' 1 in (12. 9 j) using the substitution t = 
2z/x 3 /3. If A > then 








= f^(-l)X^^, ao = 7T3 2i3-irQy u^O, (3.4) 

n=0 ^ ' 

where a n > and T(x) is the Gamma function. Hence, for the constant C we obtain 

C = a ui (l + 0(z/i)) , z/-»0. (3.5) 



To find asymptotics of 7 we use two substitutions x — v 35 and t = v sy in (12. 7p and 
(I2.8p . respectively. Then 

/■°° 2 3 ,/f _ e2 /-oo 2 3 

fl/3 + S 2 ) 2 is 



uu C / 77 —Q( s iv) ds i (3-6) 

(i/3 + S^)" 1 



OG 



where 

Q(s,z/)= / e- vt ^-^- sS Uy= I e-" n -^ t+s ?- s3 Ut. (3.7) 



Denote by L a contour in the complex s-plane which consists of the real axis with the 
interval (—1,1) replaced by the semicircle \s\ = 1, 7r ^ args is 0. Then for u < 1 



1 = uov- 1 c\ I — Q(s, 1/) rfs + 2vri Res — Q(s,i/)l. (3.8) 

\A(z/l+ S 2 ) 2 s =i,H^ + « 2 ) 2 J 

2 

One can show that Q and each of its derivatives with respect to // = U3 are uniformly 
bounded in v ^ 1, s G R. Thus, the first term in (13.81) is infinitely smooth function of /x, 
\i < 1, and 

/ " Q( a , u)ds = - f ^ + Q(i/t), i/->0. (3.9) 
7l (i/s + s 2 ) 2 Jl s 

The main term in (13. 9p can also be written in the form 

-[9M> ds = f gWHM , s + m: q, (o, 0) 

</L s </-oo s 



s 2 
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where f stands for the Cauchy principle value. The second term in (13. Sp is analytic in 
us and hence 



i- 1 
1/3 — S 



2j[i Res , 77 77 Q{s,v) = -7iiQ' s (iv3,v) 



s=th 



-TCI 



2 , 



a + bun , (3.11) 



/*oo POO 

where a = — t 2 e~3* 3 dt = —1 and b = i t(t 3 — 2)e"* 3 dt = 0. Thus, we obtain 

Jo Jo 
two formulas for computation of 7: through the principal value and as a contour integral 

Q(0,0)-Q(«,0) 



= c<rf (l + 0(i/f)V (3.12) 



where 

. 5 



c = l 9 l r m j ~2 ds - ( 3 - 13 ) 

If A < then formula (I2.12p for C~ 1 differs from (12.91) only by a sign in the exponent. 
Therefore expansion (13. 4p with (— l) n a n replaced by a n holds: 

c -^l^{0 1?-^ Mi ^t^ (3 - 14) 



Thus, asymptotics (I3.5P is valid in this case. 

The Lyapunov exponent is determined from (I2.1ip and (12.101) : 
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,1-2 /"oo -1 2 ^ /"oo o 



2u 2 J. x 0-:V ' 

foo u(^--x) poo . 

_ 2 W C / =- rfx / e -"(V-*) 

1 + X 2 

wz/ 6 / — ^ R{x,v) ax — 2uv *C I 5 , (3.15) 

J-00 (z/3 + X 2 ) 2 J -00 Z/3 + x 2 



where 



R(x,v)= e v ^ t - x ^^ t - xi Ut= e ^y-i((y+-r-^)dy. (3.16) 

</x Jo 

Denote the terms on the right hand side of (13. 15f) by I\ and I2, respectively. Then I\ is 
analogous to that in (13. 6p and gives the main contribution to 7. To determine asymptotics 
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of Ii , we use the same contour L in the complex s-plane as in (13. 8p : 

i If Z/3 — S 2 Z/3 — S 2 

Ii = ujv C < — 2 R(s, v) ds + 2ni Res — 5 ^ 

[A (z/3 + <J 2 ) 2 s=ii ,3 (z/3 + S 2 ) 2 

= wi/^C j- y rfs - 7riR' x (iu^ V) + 0(l/§ ) J 

= wz/- 1 ^ |- jf + ttz + 0(z/i) } 

_! / r° Q(0,0)-Q(s,0) ds ,1 2/ 2 \ 

= wz/ C <j4 — ^ h 0{vi)> = co-3 ^1 + 0(z/3)J ? z/ — ^ 

where c is defined by ( I3.13p . The second term I2 in (13.151) has smaller order in v 



1 f sR(s, v) ds „ z/3 — s 2 

J 2 = -2wz/~3C <^ / — ^_LJ — + 2vri Res — R(s, v) 

[Jl z/3 + s 2 s=iv h (z/3 + s 2 ) 2 

= -2u>v~*C^J ^^rfs + 7rz J R^,z/) + 0(z/l)| 

= -2o;z/-5c| -/°° - Tri Res + 7riR(iv^,v) + 0(y$] 

= -2uv^C fj°° ds + 0(z/§) j = aiO(z/i), z/ -> 0. 

Thus, asymptotics (13.121) is valid for A < as well, 
(b) Using the substitution t = 2vx we compute asymptotics of C~ l in (12. 9p 



v Jo ^ 



v 2n ~\ 



-1)« ^r(3n + |) 



7T — , V — > OO. 



n ~ n\ v 12™ ' 



Therefore, for asymptotics of C we have 
Next, we find asymptotics of (12. 8ft using repeated integration by parts: 



00. 



00 „(t+£u r d e - i[,+ - 



X 

3 



e""^)/ 2x 2(5x 2 -l) ^ s ,. 

KiT^) 1 " WTx^f + ^(1 + x 2 )* + Q(;/ } l' u ^°°- 
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This asymptotics is uniform in x. Hence, from ( 12 .7p . ( I2.8p . ( I3.2ip . and (I3.22p we obtain 
a 2 „ f°° 1-x 2 / 2x 2(5x 2 -l) ^. _ 3 .\ J 



2vu 2 (1 + x 2 ) 3 V v{l + x 2 ) 2 v 2 (l + x 2 Y 

a 2 ( I 5 



,.1-— I/" 2 + 0(0 , z/^oo. (3.23) 
ur \ lb J 

(c) The exponent in the integrand ( 12 . 1 2[) has a maximum at x = 1. Therefore, Laplace's 
method [18] immediately gives 

O- 1 = ^(1 + 0(0), ^->oo. (3.24) 
We rewrite the Lyapunov exponent as follows: 



2 3 , r°° ,/i 



(l+x 2 ) 2 

2wO / — — dx I e- u{ -- t] dt. (3.25) 

1 + X 7x 

Observe that the maximum of the exponent in the second integral is attained at t — 
l,x — —1. Therefore, Laplace's method applied to the double integral (13. 25|) (see [18] ) 
together with (I3.24p gives 

3 

Xe 2 dx e-^Ut 

-oo 1 ~l~ X J X 

= 2wC^ eT (1 + 0(^ 1 )) = w(l + Of^ 1 )), u^oo. (3.26) 

Since the integrand of the first term in (13.251) vanishes at x = —1, the same approach 
leads to the asymptotics 

/OO 1 2 3 l>CG 3 
^e v{x -~ x) dx e- u{t -- t) dt = 0{v- 2 ), u^oo. (3.27) 
oo^ + X) Jx 

Hence, in this case we have 

7 = w(l + 0(O). (3-28) 

□ 



4 Numerical results 

Since we do not estimate the remainders in asymptotic formulas (I3.ip -( 13~H]) . we compare the 
asymptotic formulas numerically with the exact ones given by (I2.7p - (l2.10p in order to determine 
the range of their validity. Figures Q] and [2] show the graphs of the normalized Lyapunov 
exponent j/lu (\2.7h as a function of v = 2u 3 /a 2 , u = \/|A[, for A > 0. As v — )■ in figure [U 
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0.02 0.04 0.06 0.08 0.1 

v 

Figure 1: Graph of the normalized Lyapunov exponent j/w of the Schrodinger equation with 
white noise potential (solid line) and its asymptotic approximation j/u — 0.3645z/~ 1 / 3 (circles) 
for small v = 2w 3 /er 2 , u = v|A[, and A > 0. 

the relative error of the asymptotic approximation of 7 becomes less than 0.7% when v < 10~ 3 . 
The relative error of asymptotic approximation ( 13. 2 j) for large values of v is less than 0.4% if 
v > 6 in figure |2j 

The graphs of the normalized Lyapunov exponent j/u H2 .101) as a function of v = 2cu 3 /a 2 , to = 
v|A| , for A < are shown in figures [3] and HI For small values of v (figure [3]) the relative error 
of the asymptotic formula (13.11) becomes less than 0.8% if v < 10~ 3 . The relative error of (13. 3p 
for large values of v (figure H]) is less than 0.7% as soon as v > 40. 
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Figure 4: Graph of the normalized Lyapunov exponent j/u of the Schrodinger equation with 
white noise potential (solid line) and its asymptotic approximation j/u = 1 (circles) for large 
v = 2u 3 /a 2 , u = a/[X[, and A < 0. 
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